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PERFORMANCE  OF  AXIAL -FLOW  TURBOMACHINES 
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SUMMARY 

A  theoretical  discussion  of  the  application  of  blade  boundary- 
layer  control  to  increase  the  efficiency  and  the  stage  pressure  ratio 
and  to  improve  off-design  performance  of  turbomachines  is  presented, 

A  method  based  on  the  potential  flow  of  a  compressible  fluid  is 
developed  for  designing  suction^  or  ejection^  slotted  blades  having  a 
prescribed  velocity  distribution  along  the  blade  and  in  the  slot.  The 
effect  of  the  boundary  layer  on  the  design  of  the  slot  and  the  effect 
of  ejecting  gas  at  staghation  pressures  and  temperatures  different 
from  free-stream  values  are  discussed. 


INTRODUCTION 

Considerable  research  has  been  conducted  on  the  application  of 
boundary- layer  control  to  isolated  airfoils.  These  investigations  prior 
to  1948  were  summarized  and  evaluated  by  Sidney  Goldstein  in  the 
Eleventh  Wright  Brothers  IVfemorial  Lecture  (reference  l) .  The  results  of 
these  investigations  suggested  certain  possible  applications  to  turbo¬ 
machine  (compressor  and  turbine)  blading,  which  were  theoretically 
investigated  at  the  NACA  Lewis  laboratory  during  the  summer  of  1950  and 
are  presented  herein. 

Boundary- layer  control  may  be  used  either  (a)  to  delay  transition 
from  a  laminar  to  a  turbulent  boundary  layer  and  thus  reduce  skin 
friction  and  heat  transfer,  or  (b)  to  prevent  or  to  delay  boundary- 
layer  separation  and  thereby  increase  the  allowable  blade  loading  and 
range  of  angles  of  attack. 

The  boundary- layer  control  to  prevent  separation  may  be  obtained 
by  either  suction  or  ejection.  The  suction  method  appears,  in  general, 
to  be  the  more  efficient,  but  the  ejection  method  would  be  necessary 
for  use  in  conjunction  with  air-cooling  of  turbine  blades. 
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Boundary- layer  suction  may  be  applied  at  a  particular  location  through 
suction  slots  or  distributed  over  an  area  by  suction  through  a  porous 
material.  Although  suction  through  porous  walls  has  certain  advantages, 
the  practical  difficulties,  such  as  clogging  and  lack  of  strength, 
make  it  appear  of  somewhat  dubious  value  for  the  blading  of  tirTbo- 
machines . 

The  applications  of  these  methods  of  boundary- layer  control  for 
increasing  the  efficiency  and  the  stage  pressure  ratio  and  improving 
the  off-design  performance  of  turbomachines  are  discussed.  A  potential- 
flow  method  of  designing  blades  for  boundary- layer  control  by  prescrib¬ 
ing  the  velocity  distribution  along  the  blade  and  in  the  slot  is 
developed  and  compared  with  the  theoretical  results  obtained  for  an 
unslotted  blade  with  a  velocity  discontinuity.  A  discussion  of  the 
effect  of  the  boundary  layer  on  the  design  of  the  slot  and  the  effect 
of  .ejecting  gas  at  stagnation  pressures  and  temperatinres  different  from 
the  free -stream  values  is  also  included. 
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SYMBOLS 

The  following  symbols  are  used  in  this  report: 

,  constants 
,M 

function  of  ^  defined  by  equation  (30) 
blade  lift  coefficient 
blade-profile  drag 
spacing  of  cascade 

complex  potential  function  (incompressible  flow) 

regular  function  of  ^ 

regular  function  of  ^ 

regular  function  of  ^ 

imaginary  part 

function  of  6  defined  by  equation  (35) 


€ 
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k 

L 

m 

n 


2q.i 

cojastant  equal  to 

constant  equal  to 

constant  defined  implicitly  ty  equations  (22)  and  (23) 
■blade-element  lift 

ratio  of  mass  flow  in  slot  to  upstream  mass  flow- 

number  determined  ly  included  trailing- edge  angle  of  Hade 


pressure 

a-uxiliary  fTUiction  of  s 

magnitude  of  dimensionless  velocity  in  compressille-flow 
plane  (ratio  of  actual  velocity  to  stagnation  velocity 
of  sound) 

R  region  in  ^-plane  defined  by  |^|>1 

Re  real  part 

s  arc  length  along  blade  surface 

u  rotor-blade-element  velocity 

v(9)  velocity  on  unit  circle  (incompressible  flow) 

w,  auxiliary  complex  variable 

z  =  X  +  iy  complex  variable  in  compressible-flow  plane  (cascade  plane) 

a  angle  of  velocity  in  compressible  flow  (measured  from 

positive  x-axis) 

p  angular  location  of  sink  on  iinit  circle 


-P 

Q(s) 

q 


r 


circulation 


0 
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p 

e 

included  trailing-edge  angle  of  blade 

C  =  1  +  in 

0 

complex  variable  in  incompressible-flow  plane 
plane) 

circle  angle  in  ^-plane 

(circle 

• 

CD 

H 

H 

W 

X 

air  angle  relative  to  axial  velocity 

p 

density 

0 

blade-row  solidity  at  mean  radius 

T 

variable  of  integration 

CP 

velocity  potential 

stream  function 

Subscripts: 

(T 

1 

upstream  of  cascade 

2 

downstream  of  cascade 

I 

entrance  to  stator  blade  row 

II 

exit  from  stator  blade  row 

III 

exit  from  next  rotor  blade  row 

a 

axial  component 

c 

compressible  flow 

i 

incompressible  flow 

m 

relative  mean  velocity 

n 

leading  edge  (nose) 

R 

relative  to  rotor 

S 

relative  to  stator  ("absolute”  values) 

• 
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t  trailing  edge  (tail) 

u  upper  surface 

The  prime  indicates  a  derivative. 


APPLICATION  OF  BOUNDARY-LAYER  CONTROL  TO  BLADING 

Although  boundary- layer  control  may  well  be  applied  to  advantage 
on  the  hub  and  the  casing  of  a  turbomachine ^  the  adequate  analysis  of 
such  an  application  involves  the  consideration  of  three-dimensional 
boundary- layer  flow  (including  the  interaction  of  wall  and  blade 
boundary  layers),  which  is  outside  the  scope  of  the  present  investi¬ 
gation  of  blade  Boundary- layer  control.  Two-dimensional  boundary- 
layer  control  on  blading  may  be  used  to  achieve  several  aims  that 
are  not  entirely  independent  but  which,  for  convenience,  are  separately 
discussed. 


Increase  of  Blade  Loading 

As  in  the  case  of  isolated  airfoils,  it  is  possible  to  increase 
the  blade  loading  considerably  by  the  use  of  boundary- layer  control. 

This  increase  results  from  the  prevention  or  the  delay  of  flow  separa¬ 
tion  and  the  resultant  stalling  of  the  blades,  and  the  greatest  gain, 
of  course,  can  be  obtained  if  the  blades  are  specially  designed  for 
the  use  of  the  boundary- layer  control  so  as  to  maintain  high  loading 
over  most  of  the  blade  section.  By  use  of  bo ixndary- layer  control,  a 
higher  blade  loading  can  be  obtained  for  a  given  inlet  Mach  number 
without  exceeding  a  given  Mach  number  on  the  blade  because  the  blade 
can  be  designed  to  have  a  uniformly  high  Mach  number  over  a  larger 
portion  of  the  upper  sirrface  of  the  blade  without  flow  separation. 

For  the  same  reason,  higher  inlet  Mach  numbers  can  be  used  with  a  given 
loading  without  exceeding  a  given  Mach  number  on  the  blade .  The  gain 
resulting  from  maintaining  a  high  velocity  over  a  larger  portion  of  the 
upper  surface  may  be  materially  reduced,  however,  by  the  increased 
velocity  over  the  lower  surface  as  a  result  of  the  practical  requirement 
of  thicker  blades  for  boundary- layer  control. 


Improvement  of  Performance  of  Later  Stage  of 
Multistage  Compressor 

The  largest  gains  from  the  increased  loading  obtainable  with 
boundary- layer  control  on  multistage  axial-flow  compressors  would  be 
in  the  later  stages  where  the  inlet  Mach  number  is  not  the  main 
limitation. 
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With  conventional  blading,  the  maximum  pressure  ratio  per  stage 
in  a  multistage  compressor  is  obtained  by  increasing  the  axial  com¬ 
ponent  of  velocity  and  maintaining  an  essentially  symmetrical  velocity 
diagram  throughout  in  order  to  maintain  the  maximum  allowable  Mach 
number  on  all  blade  elements  (reference  2) .  This  increase  in  axial 
component  of  velocity  can  only  be  obtained  by  using  a  sufficiently 
large  taper  for  the  annular  passage  to  more  than  compensate  for  the 
reduction  in  axial  velocity  due  to  the  increase  in  density.  The  use 
of  taper  large  enough  to  maintain  constant  relative  Mach  numbers, 
however,  leads  to  very  small  passages  in  the  later  stages  of  high- 
pressure-ratio  compressors  and  to  high  exit  velocities  (reference  3) . 

The  resultant  narrow  annular  passage  tends  to  produce  low  efficiencies 
in  the  later  stages  (reference  4)  and  the  high  exit  velocities  either 
produce  large  exit  losses  or  require  a  long  diffuser.  For  these 
reasons,  most  commercial  compressors  use  much  less  taper  than  required 
to  produce  constant  Mach  number  and  consequently  obtain  relatively- 
-low  Mach  numbers  and  pressure  ratios  in  the  later  stages. 

The  drop-off  in  pressure  ratio  in  the  later  stages  due  to  this 
decrease  in  relative  Ifech  number  in  the  wider  passage  could  be  pre-  p 

vented,  however,  if  the  blade  loading  could  be  Increased  in  these 
stages.  Because  of  the  lower  ^fech  numbers,  it  should  be  possible  to 
use  considerably  higher  blade  loading  without  obtaining  excessive  * 

local  Mach  numbers  on  the  blades.  With  conventional  blading,  however, 
the  blade  loading  is  limited  by  the  early  stall  of  the  blades.  The 
situation  can  be  somewhat  improved  by  the  use  of  blades  of  high  camber, 
but  the  gain  is  limited  and  the  usef\il  range  of  angles  of  attack  may 
be  reduced.  Some  further  slight  increase  in  pressure  ratio  can  be 
obtained  by  the  use  of  solidities  higher  than  conventional  limits  of 
about  1.2,  but  the  gain  is  generally  obtained  with  some  sacrifice  in 
efficiency. 

Effect  of  blade  boundary- layer  control.  -  A  much  greater  increase 
in  loading  and  stage  pressure  ratio  should  be  possible  through  the  use 
of  boundary- layer  control  on  the  rotor  and  stator  blades.  The  results 
from  isolated  airfoils  indicate  that  there  should  be  no  difficulty  in 
doubling  the  loading  obtainable  without  boundary- layer  control  with  a 
corresponding  Increase  in  press-ure  ratio.  The  effect  of  boundary- layer 
control  on  the  stage  efficiency  is  less  easy  to  evaluate  than  its  effect 
on  stage  press-ure  ratio.  The  profile  drag-lift  ratio  should  be 
decreased  because  of  the  large  increase  in  lift  coefficient  possible 
without  boundary- layer  separation  and  because  of  the  decreased  profile 
drag  resulting  from  reduced  boundary- layer  thickness  behind  the  control 
slot.  In  addition,  some  improvement  in  efficiency  might  be  expected 
from  the  fact  that  the  velocity  diagram  theoretically  most  favorable 
to  high  profile  efficiency  for  a  given  drag- lift  ratio  (symmetrical 
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diagram  vith  axial  velocity  equal  to  one-half  rotor-hlade  velocity, 
references  2  and  5)  can  he  approached  more  closely  throughout  a  multi¬ 
stage  compressor  vhen  a  high  stage  pressure  ratio  is  obtained  by  use 
of  boundary- layer  control  on  both  rotor  and  stator  blades  than  vhen 
it  is  obtained  vith  conventional  blades  by  varyiiig  the  axial  velocity 
to  give  a  constant  Mach  number  entering  all  blade  rovs.  Thus,  the 
profile  stage  efficiency  should  be  higher  than  for  conventional  blades. 
On  the  other  hand,  ”secondary-flov  losses”  may  be  increased  if  the 
empirical  relations  foimd  for  conventional  blades  (reference  6)  apply 
to  blades  vith  boundary- layer  control.  It  might  be  possible  to  reduce 
these  losses,  hovever,  by  proper  variation  of  slot  vidth  along  the 
blade  span,  by  additional  boundary- layer  control  at  critical  points, 
or  by  the  choice  of  a  suitable  relation  betveen  the  blade  spacing  and 
the  amount  of  turning  in  the  blade  rov. 

In  addition  to  the  effect  on  internal  aerodynamic  efficiency,  the 
pover  required  for  pumping  the  boundary- layer-control  air  must  be  con¬ 
sidered  in  evaluating  the  over-all  efficiency  of  the  machine.  If 
boundary- layer  control  vere  applied  to  all  blade  rovs  of  a  multistage 
compressor,  the  pover  used  in  pumping  this  air  vould  undoubtedly  have 
an  appreciable  effect  on  the  over-all  compressor  efficiency.  This 
statement  applies  even  if  the  -compressor  itself  is  used  to  pump  the 
boundary-layer  air,  as  vould  generally  be  the  case.  Because  of  this 
pumping  loss,  it  appears  desirable  to  limit  the  boundary- layer  control 
to  the  later  stages  vhere  the  largest  gains  are  possible. 

Boundary- layer  control  on  stator  blades.  -  Substantial  gains  can  be 
obtained  by  applying  the  boundary- layer  control  only  to  the  stator 
blades  of  the  latter  stages,  because  it  is  generally  necessary  or 
desirable  to  remove  the  vhirl  of  the  air  almost  completely  at  the 
discharge  from  the  compressor;  this  result  can  be  obtained  vith  a  net 
gain  in  pressure  ratio  if  the  stator-blade  loading  can  be  substantially 
increased  in  the  last  fev  stages. 

The  manner  in  vhich  this  net  gain  in  pressure  ratio  is  obtained 

can  be  seen  by  examining  typical  velocity  diagrams  at  the  mean  radius 

of  a  later  stage  vith  and  vithout  boundary- layer  control  on  stators 

(fig.  l)  .  The  subscripts  R  and  S  indicate  that  the  velocities  are 

measured  vith  respect  to  the  rotor  or  stator  reference  frame, 

respectively,  and  the  subscripts  I,  II,  and  III  indicate  the  entrance 

to  a  stator  rov,  the  exit  from  the  stator  rov,  and  the  exit  from  the 

folloving  rotor  rov,  respectively.  The  same  axial  component  of 

velocity  q^^  and  rotor-blade-element  velocity  u  are  used  in  both 

cases.  The  absolute  velocity  q^.  ^  entering  the  stator  blades  has 

l,b 

been  taken  as  the  same  in  both  cases,  but  the  air  is  shovn  turned 
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thro'ugh  a  much  larger  angle  by  the  stators  with  hovuidary- layer  control 
than  by  those  without  control^  which  results  in  absolute  velocities 
leaving  the  stators  qjj  g  that  differ  in  magnitude  and  direction. 

The  relative  velocities  entering  the  next  row  of  rotor  blades  pj 

are  therefore  also  different  and  much  higher  in  the  case  where 
boundary- layer  control  is  used  on  the  stators.  If  the  relative  Mach 
number  were  0.56  without  boundary-layer  control^  it  would  be  approxi¬ 
mately  0.75  with  control  in  this  example.  (The  txirning  in  the  stators 
should  be  limited  to  give  reasonable  relative  ^fe.ch  nvimbers  entering 
the  following  rotors.) 

The  energy  input  and  the  total  enthalpy  rise  are^  by  Euler's 
equation,  proportional  to  the  change  in  tangential  component  of  velocity 
across  the  rotor  Aq^.  The  turning  through  the  rotor  blades  in  both 

cases  of  this  example  was  determined  by  taking  Cj^  a  =  1,  where 

is  the  lift  coefficient  and  cr  the  solidity  of  the  rotor  row  of  blades. 
For  the  usual  definition  of  lift  coefficient  used  for  compressors,  this 
limitation  on  turning  gives  (reference  3) 

2|a 

q„  1 1  +  r  tan  X  „  ^ 

^a\  L 

s  neglected  and  q„  is  constant,  then 

9. 

cos 


•irI 


COS 


in,R 


C-  a  «  1  = 

J-t 


If  the  term  involving  D/l  i 


Hence,  the  total  enthalpy  rise,  with  this  limitation  on  the  turning  in 
the  rotor,  is  inversely  proportional  to  the  cosine  of  the  angle  between 
the  mean  relative  velocity  and  the  axis.  Thus,  the  use  of  boundary- 
layer  control  on  the  stators,  which  allows  this  angle  to  be  increased 
(fig.  l),  giv^s  higher  total  enthalpy  rise  and  hence  higher  total- 
pressure  ratio  than  without  the  control. 

The  process  herein  depicted  may  be  repeated  across  several  stages 
if  desired,  with  the  turning  in  the  stators  rather  arbitrarily  chosen 
to  fit  any  desired  conditions.  For  example,  the  stators  may  be 
designed  to  produce  zero  tangential  velocity  at  the  exit  of  the  last 
row  of  rotor  blades  and  thus  to  reduce  the  leaving  losses  to  a  minimum. 


Improvement  of  Turbine  Performance 

Boundaryr layer  control  appears  to  have  useful  applications  in 
conjunction  with  turbine-blade  cooling.  One  of  the  most  effective 
methods  of  cooling  the  trailing-edge  region  of  turbine  blades  is  by 
ejecting  cooling  air  at  or  near  the  trailing  edge  of  the  blade.  This 
ejected  air  can  be  used  at  the  same  time  as  an  effective  method  of 
boundary- layer  control  for  increasing  the  blade  loading  and  thereby 
reducing  the  total  blade  area  for  a  given  power  output.  The  reduction 
in  blade-surface  area  may,  in  t^Irn,  be  expected  to  reduce  the  amount 
of  cooling  required  to  maintain  a  given  blade  temperature .  A  some¬ 
what  different  application  of  boundary- layer  control  would  be  made  for 
turbines  requiring  very  hi^  work  per  potind  of  gas  with  limited  rotor- 
blade  speed,  where  negative  reaction  may  be  required.  The  flow  throxigh 
the  tiirbine  rotor  is  then  similar  to  that  throu^  a  typical  compressor 
rotor  and,  if  the  blade  loading  is  to  be  high,  boundary- layer  control 
may  be  required  to  prevent  stalling  of  the  blades.  In  this  case,  as 
in  the  more  usual  case,  boundary- layer  control  by  ejection  can 
effectively  be  used  in  conjunction  with  blade  cooling. 


Improvement  af  Off -Design  Performance  of  Turbomachines 

In  a  multistage  turbomachine  designed  for  a  high  pressvire  ratio, 
the  performance  may  be  very  poor  when  operating  at  equivalent  speeds 
and  weight  flows  far  removed  from  the  design  values.  The  loss  in 
efficiency  is  particularly  serious  in  multistage  axial-flow  compressors 
when  operating  at  below  design  speed,  and  is  a  result  of  extreme  mis¬ 
matching  of  the  different  stages.  The  first  stages  must  operate  at 
high  angles  of  attack  and  may  stall,  whereas  the  last  stages  may  be 
operating  as  a  turbine  (which  reduces  the  over-all  pressure  obtained) 
or  may  even  be  subject  to  a  negative  stall  of  the  blades. 

Boundary- layer  control  may  be  used  to  improve  this  situation  in 
several  ways.  In  the  first  place,  airfoils  properly  designed  for 
boundary— layer  control  are  capable  of  withstanding  a  much  wider  range 
of  angles  of  attack  than  those  without  control  (reference  7). 

Another  method  of  improving  the  off-design  performance  with 
boundary- layer  control  consists  in  varying  the  turning  angle  and  the 
circulation  for  a  given  angle  of  attack  by  varying  the  amount  of  air 
sucked  off  for  control.  In  general,  increasing  the  suction  throvigh  a 
slot  near  the  tail  edge  on  the  upper  surface  increases  the  circulation 
(reference  8) .  This  variation  of  the  suction  quantity  can  be  used 
very  effectively  in  the  later  stages  of  a  multistage  compressor  where 
boundary- layer  suction  has  already  been  suggested  for  improving  the 
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design  performance.  Similar  effects  are  obtained  ty  varying  ejection 
air  in  the  ejection-control  method  for  perfect  fluid,  but  the  analysis 
of  the  effect  for  actual  fluid  is  not  as  clear. 

The  tendency  of  the  later  stages  to  perform  little  useful  work  at 
below-design  speed  due  to  a  large  reduction  in  angle  of  attack  can  be 
counteracted  by  an  increase  in  the  circulation  by  using  increased 
suction  at  the  reduced  speeds.  This  benefit  is  obtained  even  when  the 
suction  is  applied  only  to  the  stator  blades  because  an  increased 
turning  through  the  stator  blades  results  in  an  increased  angle  of 
attack  on  the  following  row  of  rotor  blades  and  hence  an  increase  in 
energy  addition  through  the  rotor  blades. 

Boundary- layer  control  has  not  been  particularly  recommended  for 
improving  the  design  performance  of  the  first  stages  of  axial-flow 
compressors  because  of  the  high  relative  entrance  Mach  numbers  required 
on  these  stages  to  obtain  high  mass  flow  and  high  rotor  speed.  The 
high  Mach  numbers  require  thin  blades  with  relatively  low  loading,  in 
which  case  the  gain  from  boundary- layer  control  is  small.  It  might  be 
desirable,  however,  to  employ  boundary- layer  control  on  these  stages 
during  the  low-speed  operation  of  the  compressor  in  order  to  postpone 
the  stalling  of  these  blades  at  the  high  angles  of  attack  at  which  they 
will  then  operate.  Suction  near  the  leading  edge  is  fairly  effective 
in  postponing  the  stalling  of  thin  airfoils  (references  1  and  9  to  11) 
and  might  therefore  be  used  to  advantage  in  improving  the  performance 
of  the  first  stages  at  low  compressor  speeds. 


THEORY  OF  BLADE-DESIGH  METHOD  WITH  SLOT 

In  reference  12,  Lin  shows  that,  if  the  pressure-density  relation 
is 


P  = 


(1) 


then  the  compressible  potential  flow  about  a  cascade  of  blades  can  be 
obtained  by  transforming  an  incompressible  flow  about  the  vinit  circle 
with  the  unit  circle  mapping  into  the  cascade  blades.  This  transforma¬ 
tion  is  extended  to  give  the  compressible  potential  flow  past  a  cascade 
of  blades  with  suction  or  ejection  slots  in  the  blades  by  Introducing  a 
sink  or  source  on  the  unit  circle  and  a  pole  in  the  derivative  of  the 
mapping  function  at  the  same  point. 

The  complex  potential  fxxnction  F(^)  for  the  incompressible  flow 
about  the  unit  circle  due  to  complex  sources  at  ^  =  ±e  (k  real 


NACA  031  2371 


11 


and  positive)  and  a  sink  or  source  on  the  circle  at  ^  =  e^^  (p  real) 
is 

F(0  =  A  loggC^  “  +  A  logg(^  -  e"^)  + 

B  loSe(^  +  e^)  +  B  logg(C  +  e"^)  +  2  M  logg(^  -  e^^)  +  C  (2) 

where  the  har  indicates  the  complex  conjugate^  M  is  a  real  constant 
(positive  for  a  source  and  negative  for  a  siiik:)  and  A,  and  C  are 
complex  constants  with 

Re  A  =  (-  Re  B  -  M)  >  0  (3) 

The  mapping  "between  the  z-plane  and  the  ^-plane  defined  hy 


gives  a  compressible  flow  with  a  linear  pres sure- volimie  relation  past  a 
cascade  of  slotted  blades  with  the  velocity  potential  cp^  and  stream 
function  in  the  z-plane.  given  by 

cpg  +  i>lfg  =  F(0  (5) 


provided  that  g(^)  is  chosen  to  satisfy  the  following  requirements 

(a)  g(^)  is  regular  in  region  R  defined  hy  |  >1 

(h)  g(^)^  0  in  R  except  where  F*.(^)  =  0  (the  order  of 
the  zero  not  to  exceed  1) 


dz  =  0  along  a  path  enclosing  the  unit  circle  and 

h 

excluding  the  points  ^  =  ±e 
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The  magnitude  q,  and  direction  a  of  the  dimensionless 
velocity  at  any  point  in  the  z-plane  are  given  by 


2q 


la 


2k, 

e  ) 


1  +  q' 


(■■•r) 


(7) 


The  proof  that  the  differential  equations  for  the  velocity  poten¬ 
tial  and  stream  function  are  satisfied  and  that  the  required  boiindary 
conditions  are  obtained  in  the  z-plane  by  use  of  this  transformation  is 
the  same  as  the  proof  given  in  reference  12.  In  the  physical  interpre¬ 
tation^  part  of  the  suction  slot  is  considered  as  being  ducted  out  along 
the  span  of  the  blade. 


The  preceding  transformation  can  be  extended  to  give  additional 
suction,  or  ejection,  slots  on  the  blade  by  adding  more  sinks  or 
sources  on  the  xmit  circle  and  by  putting  poles  in  the  derivative  of 
the  mapping  function  at  the  same  points.  In  order  to  use  this  trans¬ 
formation  in  designing  a  cascade  of  blades  with  suction  or  ejection 
slots  and  a  prescribed  velocity  distribution  along  the  blade,  the 
prescribed  conditions  are  used  to  select  a  suitable  Incompressible  flow 
about  the  unit  circle  and  to  determine  the  function  g(0  ^  manner 

similar  to  that  used  for  a  cascade  of  unslotted  blades  (references  13 
and  14). 

Tlie  prescribed  conditions  are  the  velocity  distribution  on  the 

blade,  the  upstream  velocity  q^  the  downstream  velocity  q^ 

and  the  ratio  m  of  mass  flov;  in  the  slot  to  inlet  mass  flow  (m  is 
negative  for  auction,  positive  for  ejection).  Because  the  prescribed 
velocity  does  not  determine  the  circulation,  either  the  circulation  or 
the  location  of  two  equipotential  points  in  the  slot  must  also  be 
assigned.  The  upstream  and  the  dowistream  velocities  are  related  by  the 
continuity  equation 


^2 


cos  ag 


(l4m)  cos 


(8) 


Because  the  suction  case  and  the  ejection  case  differ  only  in  the 
signs  of  M  and  m  and  in  the  location  of  the  upper- surface  stagna¬ 
tion  point  relative  to  the  slot,  only  the  suction  case  is  considered 
in  detail.  For  convenience,  the  axis  of  the  cascade  is  taken  along 
the  y-a;ds  and  the  flow  is  from  right  to  left  (fig.  2). 


2116 


MCA  !M  2371 


13 
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■» 


Flow  in  Circle  Plane 

The  flow  of  an  incompressible  fluid  about  the  unit  circle  is 
selected  by  determining  the  constants  k,  P,  and  k  in  the  complex 
potential  function  F(^)  from  the  given  conditions.  The  constants  A 
and  B  are  determined  from  the  upstream  and  the  downstream  velocities 
and  the  circulation}  and  p  and  k  are  determined  by  the  range  of 
potential  between  stagnation  points. 

Potential  and  circulation  on  blade.  -  The  magnitude  of  the 
dimensionless  velocity  q  along  the  surface  of  the  blade  and  the 
slot  is  given  as  a  function  of  the  arc  length  q  =  Q(s)^  where  the 
total  arc  length  is  infinite  and  the  trailing-edge  stagnation  point 
s^  has  been  taken  at  zero  (fig.  3) 

Q(s)  =  -  q(s) 

Q(s)  =  q(s) 

Q(s)  =  -  q(s) 

Q(s)  =  q(s) 

where  and  are  the  upper  surface  and  the  nose  stagnation 

points,  respectively,  then  the  velocity  potential  cp^(s)  is  given  by 

s 

cp^(s)  a  r  Q(s)  ds  (10) 

"^0 


If  Q(s)  is  defined  hy 


-oo<  s<  s.^ 
s^<  s<  0 
0<  s<  Sj^ 
Sn<  s< +  00 


MS) 


The  circulation  cannot  be  determined  from  the  prescribed  velocity 
because  the  integral 


+  00 

J'  Q(s)  ds 


is  an  indeterminate,  improper  integral.  Hence,  either  the  circulation 
or  the  location  (as  a  function  of  s)  in  the  slot  of  two  points  a  and 
b  having  the  same  potential  must  be  assigned  (figs.  2  and  3).  When 
the  location  of  points  a  and  b '  is  assigned, 

^  J' 


(11) 
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With  the  circulation  known;  the  spacing  d  of  the  cascade  is 
given  by 


d  = 


sin  -  qg  sin  a2 


(12) 


Determination  of  A;  M;  and  B.  -  The  value  of  d  from  equa¬ 
tion  (12)  is  used  to  evaluate  A  and  B  because  the  spacing  is  also 


given  by  the  value  of 
that  iS; 


taken  along  a  path  around 


k 


e 


or 


(13) 


Substituting  dz  from  equation  (4)  in  equation  (13)  and 
integrating  give 


f 


H 


id  =  2jti 


=  2ni 


6(23 


(l  -  4  i(e3 


A  ia_  AKt  e 

—  e  ^  +  — — 
Ki  4 


la- 


(14) 


where 


2q. 


K,  =  - , 

1  H-a/l  + 


The  bracketed  expression  must  be  a  real  number,  hence, 

2 


4+K- 


Im  A  = 


4-K 


■2  Re  A  tan 


1 


(15) 


or 


Im  A  =  +  Re  A  tan 


(16) 
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# 


* 


% 


Substitution  of 


Im  A  from  equation  (l5) 


Re  A  = 


4K^  cos  a. 


4  + 


in  equation  (14)  gives 

2jt 


(17) 


or 


Re 


q^^  I  cos 

's/iT^ 


d 

2jt 


Hence,  A  is  completely  determined  by  equations  (15)  and  (17)  or 
(16)  and  (18). 


(18) 


The  ratio  m  (which  is  negative  in  the  suction  case)  of  the  mass 
flow  removed  to  total  flow  is  given  and  consequently  e 

M  i=!  m  Re  A  (l9) 

Equation  (3)  may  be  written 

Re  B  =  -  Re  A  -  M  (20) 

and  Re  B  is  thus  obtained  from  the  preceding  values  of  Re  A  and  M. 
The  value  of  Im  B  is  determined  from  the  spacing  by 


where 
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Therefore 


4+K< 


Im  B  = 


4-K, 


Re  B  tan 


^-4. 


1  + 


Re  B  tan  a. 


qg  - - 2 

and  B  is  completely  known  from  equations  (20)  and  (21)  . 


(21) 


Determination  of  k  and  p.  -  After  A,  and  B  are  known,  k 
and  B  are  determined  by  the  condition  that  the  maximums  and  the  mini- 
mirniF!  of  the  potential  on  the  circle  must  equal  those  of  the  potential 
on  the  blade;  that  is, 


%(sj  -  %(0)  =  F(e''®^)  -  F(e''®^) 

(22) 

(23) 

where  0^,  0^,  and  0^^  are  the  upper  surface,  the  trailing  edge,  and 

the  nose  stagnation  points,  respectively.  This  system  of  equations  in 
k  B  always  has  solutions,  but  the  equations  cannot  be  solved 

analytically  and  some  numerical  method  must  be  used.  One  method  is  to 
obtain  a  plot  of  k  as  a  function  of  B  values  that  satisfy  equa¬ 

tion  (22)  by  assigning  a  few  values  to  B  and  determining  the  corre¬ 
sponding  values  of  k.  Similarly,  obtain  a  plot  of  k  as  a  function 
of  B  fon  values  that  satisfy  equation  (23) .  The  intersection  of  these 
curves  gives  the  desired  values  of  k  and  B* 

With  k  B  determined,  the  flow  about  the  circle  is  known. 

The  potential  cpp(0)  on  the  circle  is 

cPi(0)  =  Re  =  -  2  Re  A  tanh"!  + 


(Un  A  +  Im  B)  tan'l  +  (Im  A  -  Im  B)  tan’l  t 


tanh  k 

M  1  -  cos  (0  -  B) 

^  ®e  cosh  k  +  cos  0 


+  Ra  C 


(24) 
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and  Re  C  is  chosen  to  mahe  =  0.  The  angle  convention  is 


- 


s  in  0  ^  jt 
sinh  k  2 


and 


tan"^  -I—,  — is  taken  in  the  same  quadrant  and  direction  as 
tanh  k 

The  velocity  on  the  circle  v(0)  is 


e. 


v(e)  = 


2  Re  A  sin  0  cosh  k  + 


cosh  2k  -  cos  20 

(Tth  a  -  Im  B)  cos  0  sinh  k  +  (im  A  +  Im  B)  sinh  k  cosh  k 

,  0-p  ,  sin  0 

V  ^  2  cosh  k  +  cos  0 


(25) 


% 


Function  g(^) 

The  function  g(0  can  he  computed  for  points  on  the  unit  circle 
by  using  the  prescribed  velocity  on  the  blade  and  the  velocity  on  the 
circle  to  determine  the  real  part  of  g(^)  .  The  imaginary  part  could 
then  be  computed  by  Poisson’s  integral.  Because  of  the  restriction 
imposed  by  the  prescribed  conditions^  g(^)  is  obtained  in  a  slightly 
different  manner ^  as  shovn  in  the  following  sections. 


Blade  with  pointed  trailing  edge.  -  If  the  blade  is  to  have  a 
pointed  trailing  edge^  then  g^"^  must  vanish  at  the  trailing-edge 


stagnation  point  on  the  circle,  ^  =  e 
as 


i0. 


Hence,  g(0  ^^7  written 


s(S)  = 


where  f(^)  is  regular  on  and  outside  the  unit  circle  and 


(26) 


(27) 


1. 


where  e  is  the  included  trailing-edge  angle  of  the  hlade. 
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Values  of  g(±e^) .  -  Because  the  velocities  upstream  and  down¬ 
stream  are  prescribed,  the  values  of  g(0  determined 

from  equation  (7); 


g(e^) 


k-)  ^  (e^  -  e^P) 


K- 


M28) 


g(.ek)  _  (el^  ^  ) 


In  order  that  g(^)  have  these  values,  f(^)  is  written  in  the  form 


f(0  =  C(^)  + 


(e2k  _  ^2)  (^2  _  e~2k) 


H(0 


(29) 


where 


k 


5  11  - 


e 

T 


-2B(e^  +  e^P)e^°’^ 
Kp(l  +  e^^t-^)"^ 


(30) 


and  H(^)  is  regular  on  and  outside  the  unit  circle  with 

lim  ^  H(^)  =  0 
>co 


(31) 


This  restriction  on  H(C)  is  imposed  by  the  regularity  of  f (^)  .  By 
use  of  equations  (26)  and  (29) , 


g(0 


le. 


n 


®  e[c(?)  +  (e21^-C^)(^2_e-2k)^-2 


(32) 


and  g(0  will  be  known  when  H(U  is  determined.  For  the  computation 
of  the  blade  shape,  only  the  values  of  g(0  on  circle  are  needed. 
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Determination  of  Re  H(0  on  circle.  «  By  equation  (5),  the 
potentials  cp^(s)  and  C(:^(0)  are  equal  at  corresponding  points.  Thus, 

hy  matching  these  potentials,  a  correspondence  is  established  between 
points  along  the  blade  arc  length  and  the  circle  angle:  s  =  s(0).  By 
use  of  this  correspondence,  the  magnitude  of  the  velocity  prescribed 
along  the  blade  is  obtained  as  a  function  of  the  circle  angle: 

=  <l(^)  •  Hence,  from  equation  (7), 


2  q(e) 

F'(e^®)(ei20  . 

■  .21)  ( 

1  _  e^(P"®))l 

1  +a) 

^1  +  q^(0) 

|g(ei0)' 

1 

Substituting  g(^)  from  equation  (32)  with  ^  =  e^^  and  replacing 
jF^(e^^)j  by  the  velocity  v(0)  (equation  (25))  in  equation  (33)  give 


> 


2q(9) 


cos 


e^|v(0)|(2  cosh  2k-2  cos  20)^^^  2  sin 

-|n/2  [rc  C(ei9)  +(2  cosh  2k  -  2  cos  20)Re  H(ei0)3 
( 0-(;-0)J  e 


or 


v(0)t(2  cosh  2k -2  cos  20)^/^  2 


log. 


.  p-e 

sin 


Re  E(e^®)  = 


K{0)j^2  -  2  cos  (0^-0^ 


n/2 


Re  C(e^®)  +  k 


2  cosh  2k  -  2  cos  20 


(34) 


where 


K(0)  =  - . -  (35) 

1  +a/i  +  q^(e) 

j  Restrictions  on  Re  H(e^^) .  -  By  equation  (31) , 

lim  ^  H(0  =  0 

•  5-^0 
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hence  H(0  must  have  the  form 


(vhere  hg,  h^^  .  .  .  are  complex  constants)  vhich  for  points  on  the 
unit  circle  becomes 


H(e^®)  =  P.e  +  i  Im  H(e^®)  ^ 


i0 


10' 


,120 


,130 


0 


=  E  (Re  h.  cos  j0  +  Im  h,:  sin  j0)  + 
3=2 


1  E  (im  hi  cos  j0  -  Re  hi  sin  j0)  (36) 

J  =  2  J  J 

Consequently, 


/ 

Re  H(e^®) 

d0  =  0 

(37) 

0 

2jt 

Re  H(e^®) 

f 

cos  0  d0  =  0 

(38) 

0 

2rt 

Re  H(e^®) 

/ 

sin  0  d0  =  0 

(39) 

0 


The  values  of  Re  H(e^®)  from  equation  (34)  must  satisfy  equa¬ 
tions  (37)  to  (39) .  If  these  equations  are  not  satisfied,  the  values 

of  Re  H(e  )  must  he  adjusted  until  the  conditions  are  satisfied, 
which  modifies  the  prescribed  velocity,  tfethods  for  adjusting 
i0 

Re  H(e  .)  are  given  In  reference  14.  These  adjustments  must  he  made 

In  such  a  manner  that  Re  H(e^®)  still  satisfies  the  restriction 
Imposed  by  equation  (6d),  which  may  be  written  In  the  form 


* 


¥ 


¥ 
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5 


y 


% 


gr(2  cosh  2k  -  2  cos  26)  Re  H(e^®) 


e^ 

v(e) 

(2  cosh  2k  -  2  cos  26)^/^ 

1  .  P-0 

2  -  2  cos  (6^-6) 

n/2  gR6  c(e 

.16) 

(40) 


Detennimtion  of  Im  H(e^^) .  -  After  obtaining  Re  H(e^®),  which 
satisfies  equations  (37)  to  (40),  the  fimctlon  Im  H(e^®)  is  given  by 
Poisson's  integral 


Im  H(e^®) 


^2n 


dT 


(41) 


where  the  constant  term  in  the  integral  has  been  taken  as  zero  so  that 


as  required  by  equation  (36) .  Hence,  H(e^®)  is  now  determined  by 
E(ei®)  =  Re  H(e^®)  +  i  Im  H(e^®) 

and  substitution  of  these  values  of  H(e^®)  in  equation  (32)  gives  the 
values  of  g(^)  for  points  on  the  circle  ^  =  e^®. 


Blade  Coordinates 


From  the  preceding  values  of  g(e^®),  the  blade  coordinates  are 
obtained  by  integration  around  the  circle  of  equation  (4)}  that  is. 
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which,  on  replacing  F'(e^^)  hy  v(0)  and  writing 


g(e^®)  =  g3^(0)  e 


ig2(®) 


reduces  to 


z  = 


g^(e)  Cl  -  (e2i0  .  . 


2 

^  g3_(e) 


v{e)  ^  _  ^i(p.eg  ^^2ie  _  ^2k, 


ILLUSTRATIVE  EXAMPLE 


i(e  +  J +  g2) 


d0 

(42) 


In  order  to  compare  an  example  of  a  slotted  blade  with  a  previously 
computed  unslotted,  Griffith  type  blade,  which  had  an  abrupt  decrease 
in  the  upper  surface  velocity  (figs.  4  and  5),  the  same  upstream 
velocity,  turning  angle,  spacing,  and  circulation  were  chosen  for  this 
example : 


=  0.68215 
a,^  =  135° 

ag  =  180° 
r  =  0.5058 
d  =  1.0488 

and  a  similar  velocity  distribution  was  prescribed  along  the  blade. 
The  amount  of  suction  was  taken  to  bo  1  percent  of  the  mass  flow 
(m  =  -O.Ol).  The  integrals  in  equations  (36)  to  (40)  were  evalmted 
by  the  methods  given  in  reference  14.  In  order  to  obtain  sufficient 
accuracy  in  the  nei^borhood  of  the  slot,  when  the  coefficients  given 

i0 

in  reference  (14)  are  used.  Re  H(e  )  was  written  in  the  form 
Re  H(ei®)  =  Re  H^(e^®)  +  Re  HgCe^®) 
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where 

Re  %(e^®) 

Re  H^(ei®) 

Re  HgCe^®)  =  Re  H(e^®)  for  p  -  20°  <  0  <  p  +  20° 

Re  H„(e^®)  =  0  for  all  other  values  of  0 

u 

The  integrals  involving  Re  H-,  (e^Q)  were  evaluated  using  80  equally 
spaced  values  of  0  and  the  integrals  involving  Re  H  (e  )  were 

obtained  hy  conformally  mapping  the  imit  circle  onto  a  unit  circle  hy 

w  =  ^  -  0.95  e^^ 

0.95  e"^^  ^  -  1 

and  using  80  equally  spaced  values  of  the  new  circle  angle. 

Adjustments  to  Re  H(e^^)  changed  the  prescribed  velocity  somewhat 
and  the  final  velocity  distribution  is  shown  in  figure  6.  The  computed 
blade  shape  is  given  in  figure  7 .  The  trailing  edge  of  the  slot  has  a 
very  small  radius  of  curvature  and;  in  order  to  obtain  a  more  practical 
blade;  the  velocity  was  represcribed  in  this  neighborhood  (fig.  8)  and 
also  a  lower  velocity  was  prescribed  in  the  slot.  The  resulting 
velocity  distribution  and  blade  are  shown  in  figures  9  and  10; 
respectively. 

DISCUSSION  OF  SLOT  DESIGN  METHOD 

The  design  method  with  a  slot  directly  applies  to  potential  com¬ 
pressible  flow  with  a  linear  pressure-volume  relation.  The  inter¬ 
pretation  given  in  reference  14,  which  consists  in  multiplying  the 
velocities  for  a  linearized  pres sure- volume  relation  by  a  constant 
factor;  can  be  used  to  obtain  an  approximate  solution  for  a  gas  obeying 
the  usual  isentropic  relations  and  satisfying  irrotationality  and  con¬ 
tinuity  between  upstream  and  downstream  flow.  The  limitation  to  poten¬ 
tial  floW;  however;  is  a  more  serious  restriction  because  most  of  the 
slot  flow  is  boundary- layer  flow  with  total  pressures  and  temperatures 


=  0  for  P  -  20°  <  0  <  P  +  20° 

=  Re  H(e^®)  for  all  other  values  of  0 
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differing  from  those  of  the  main  flow.  Thus,  the  method  cannot  he 
directly  applied  to  actual  flows.  A  reasonable  method  of  allowing  for 
the  bo\andary  layer  in  the  case  of  suction  would  be  to  displace  the 
upstream  boundary  of  the  slot  by  a  little  more  than  the  displacement 
thickness  of  the  boundary  layer  in  order  to  allow  the  slot  to  pass  the 
desired  mass  flow. 

The  modification  of  the  potential-flow  solution  to  represent 
actual  ejection  control  is  more  complicated.  Even  without  the  effect 
of  viscosity,  the  ejection  process  may  differ  from  the  potential-flow 
solution  because  the  total  pressure  and  temperatmre  of  the  ejected 
gas  need  not  be  the  same  as  that  of  the  main  flow.  Then  a  discon¬ 
tinuity  in  both  temperature  and  velocity  (but  not  pressure)  would 
generally  be  produced  where  the  two  flows  come  together.  'When 
viscosity  is  considered,  the  discontinuity  surface  disappears  and  will 
be  replaced  by  a  wake  region  where  the  boundary  layers  from  the  blade 
and  the  upstream  side  of  the  slot  merge. 

The  nonviscous  case  can  be  analyzed  approximately  on  a  one- 
dimensional  basis  by  applying  Berno-ulli's  equation  to  stream  tubes  on 
either  side  of  the  surface  of  discontinuity  and  making  use  of  the 
continuity  of  pressure  across  the  surface  of  discontinuity.  The  viscous 
case  is  considerably  more  difficult  and  must  be  treated  by  boundary- 
layer  techniques  that  are  beyond  the  scope  of  the  present  investigation. 
The  general  approach  to  obtaining  the  actual  flow  with  ejection  would 
be,  however,  to  start  with  the  potential-flow  solution  obtained  in 
the  present  report  and  then  to  modify  it  in  succession  to  take  account 
of  different  total  press\ares  and  temperatures  and  then  modify  that 
solution  to  account  for  viscous  and  heat- transfer  effects. 

For  low  suction  or  ejection  flows  (l  percent  or  less  of  upstream 
flow),  the  blade  shape  obtained  by  the  method  of  reference  14  using  a 
discontinuity  in  the  velocity  is  essentially  the  same  as  the  blade 
shape  obtained  by  the  method  of  this  report  except  in  the  vicinity  of 
the  slot.  Consequently,  considerable  time  could  be  saved  in  designing 
blades  with  low  suction  or  ejection  flows  by  using  the  method  of  refer¬ 
ence  14  to  obtain  the  general  blade  shape  and  then  fairing  into  the 
blade  a  previously  calculated  (or  known)  slot  that  is  capable  of 
handling  the  flow. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  December  13,  1950. 
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Fioure  7.  -  First  slotted  blade  with  suction.  Ratio  of  mass  flows,  m,  7O.OI; 
cascade  spacing,  d,  1.0488;  upstream  of  cascade:  magnitude  of  dimensionless 
velocity,  Qj ,  0.68215;  angle  of  velocity,  a|,  135®;  downstream  of  cascade: 
magnitude  of  dimensionless  velocity,  Qg.  0.42930;  angle  of  velocity,  a^. 
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Figure  10.  -  Second  slotted  blade  with  suction.  Ratio  of  nass  flows,  tn,  -0.01;  cascade 
spacing,  d,  1.0488;  upstream  of  cascade:  magnitude  of  dimensionless  velocity,  Q] ,  4 

0.68215:  angle  of  velocity,  «|,  135®;  downstream  of  cascade:  magnitude  of  dimension¬ 
less  velocity,  03,  0.42930;  angle  of  velocity,  a2,  180®. 
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